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ABSTRACT 

Most of the baryons from galaxies have been "missing" and several studies 
have attempted to map the circumgalactic medium (CGM) of galaxies in their 
quest. We report on X-ray observations made with the Chandra X-ray Observa- 
tory probing the warm-hot phase of the CGM of our Milky Way at about 10 6 K. 
We detect O VII and O vm absorption lines at z = in extragalactic sight lines 
and measure accurate column densities using both Ka and K/3 lines of O VII. We 
then combine these measurements with the emission measure of the Galactic halo 
from literature to derive the density and the pathlength of the CGM. We show 
that the warm-hot phase of the CGM is massive, extending over a large region 
around the Milky Way, with a radius of over 100 kpc. The mass content of this 
phase is over ten billion solar masses, many times more than that in cooler gas 
phases and comparable to the total baryonic mass in the disk of the Galaxy. The 
missing mass of the Galaxy appears to be in this warm-hot gas phase. 
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1. Introduction 

We have known for a while that the baryonic mass of galaxies, including that of our 
own Milky way, falls short of what is expected for their total mass ((Sommer-Larsen 2006; 
Bregman 2007) and references therein). This "missing" mass is either ejected from galaxies 
into the intergalactic medium (IGM) or still resides in the circumgalactic medium (CGM). 
Our galaxy is a member of the Local Group of galaxies, and it is possible that the matter 
ejected from the Galaxy resides in the Local Group. Either in the CGM or in the Local 
Group medium (LGM), the baryonic mass is expected in a warm-hot gas phase at temper- 
atures between 10 5 — 10 7 K. This is probably the reason why earlier studies probing only 
cooler gas could not account for it. The warm component of this phase at 10 5 — 10 6 K 
has been observed in the UV traced by absorption lines of ionized metals, in particu- 
lar Ovi (Nicastro et al. 2003; Sembach et al. 2003). Recent observations with the Hubble 
Space Telescope (HST) have shown that the CGM around star-forming galaxies is large (150 
kpc) and the mass of the warm phase traced by O VI exceeds that of the gas in the galax- 
ies themselves (Tumlinson et al. 2011; Tripp et al. 2011). Using other absorption lines, e.g. 
Si IV and C IV in the UV (Lehner & Howk 2011) found that there is a large reservoir of warm 
ionized gas around our Galaxy as well. 

The hotter phase of the warm- hot gas, at temperatures 10 6 — 10 7 K can be probed by 
even more highly ionized metals. The dominant transitions from such ions lie in the soft- 
X-ray band and indeed, absorption lines due to O VII and O vm at redshift zero have been 
detected toward extragalactic sight-lines by the Chandra X-ray Observatory and the XMM- 
Newton Observatory (Nicastro et al. 2005; Williams et al. 2005; Williams et al. 2007; Bregman et al. 200' 
The distribution, spacial extent, and mass of the warm- hot gas provide important constraints 
to the models of large scale structure formations (e.g., (Cen & Ostriker 1999)) and/or out- 
flows/inflows from galaxies (e.g. (Stinson et al. 2011)), but these parameters are difficult to 
measure in part due to the limited spectral resolution of Chandra and XMM-Newton and in 
part because of the inherent difficulty in using the absorption line studies alone. Here we use 
both absorption and emission observations to derive the physical properties of the warm-hot 
plasma. We show that the X-ray observations probe million degree gas, with low density, 
extending over 100 kpc and having mass over ten billion solar masses. This is several times 
more than previously found in the CGM of the Milky Way. Alternatively, the warm-hot gas 
we probe is from the extended LG medium. 
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2. Sample selection & Data reduction 

We search for and measure the absorption lines from highly ionized gas at z ~ 0, as 
seen in the spectra of most of the extragalactic sight-lines toward bright active galactic 
nuclei (AGNs). We use high signal-to-noise ratio observations with Chandra High Energy 
Transmission Grating (HETG) or Low Energy Transmission Grating (LETG) with focus on 
O VII and O vm absorption lines at A = 21.602 A and A = 18.967 A respectively. 

We began by selecting all the Chandra grating observations of AGNs that were publicly 
available as of 2011 October 30 and that had exposure time of at least 100 ks. This resulted in 
multiple observations of 50 targets with both grating spectrometers aboard Chandra (HETG 
and LETG). 

We reduced the data using the standard Chandra Interactive Analysis of Observa- 
tions (CIAO) software (v4.3) and Chandra Calibration Database (CALDB, v4.4.2) and fol- 
lowed the standard Chandra data reduction threads. For the Chandra ACIS/HETG and 
ACIS/LETG observations, we co- added the negative and positive first-order spectra and 
built the effective area files (ARFs) for each observations using the fullgarf CIAO script. 
Those pertaining to the ACIS/LETG observations were corrected for the ACIS quantum 
efficiency degradation. Unlike ACIS, the HRC does not have the energy resolution to sort 
individual orders, and each spectrum contains contributions from all the diffraction orders. 
For the HRC/LETG observations we used the standard ARF files for orders 1 through 6 and 
convolved it with the relevant standard redistribution matrix file (RMF). For the targets 
with multiple observations, we added the grating spectra and averaged the associated ARFs 
using the ciao script add_grating spectra, to increase the signal to noise ratio of the spectra. 
We only add the observations with the same instrumental configuration and the observations 
made with different instruments are analyzed separately. 

Of the 50 sight lines initially selected, 29 have good enough signal-to-noise ratios (S/N) 
near 21.602 A to detect Ovu absorption lines. The strong Ovil Ka absorption lines near 
21.602 A are clearly visible in 21 out of 29 sources with good signal-to-noise spectra; thus 
the covering fraction of the O VII systems is 21/29 = 0.72. In 30% of the sources significant 
O vm Ka lines near 18.967 A are also observed. Here we consider a sub-sample of 8 targets 
(Table 1), where both Ovil and Ovill Ka local absorption has been confidently detected 
(we have not included 3C273 in this sample of 29 sources because the z = absorption 
lines in this sight-line might be from a nearby supernova remnant). The detailed analysis 
of the complete sample will be presented in a forthcoming paper (Gupta et al. 2012, in 
preparation). The local (z ~ 0) Ovil and Ovm absorption lines in 3 of our 8 targets have 
been reported previously by other authors, but to ensure the consistency of data analysis, we 
reanalyzed all the data and obtained the fit results independently; for the other five targets 



-4- 



we present new detections of z ~ lines. 

3. Analysis 

3.1. Continuum, Intrinsic Absorption, and Emission 

The HETG and LETG spectra are binned to 0.01 A and 0.025 A respectively and 
analyzed using the CIAO fitting package Sherpa. Since we are interested in oxygen absorption 
lines we fit the spectral continuum in the 17 — 23 A range with a simple power law and a 
Galactic absorbing column density (N#, from (Dickey et al. 1990)). For the targets Mrk290 
and NGC3783 we need to add a black-body component to account for the soft excess, that 
could come from a standard optically-thick accretion disk. In our sample, the photon index 
of the power-law varies from 1.4 — 2.1 with average value of 2.0. 

Our targets are nearby Type 1 AGNs, which have their own intrinsic absorption and 
emission features. We carefully study all the spectra to confirm that none of the AGN 
intrinsic features contaminates the local (z = 0) absorption lines. Except for one source, 
NGC4051, the intrinsic lines are sufficiently red-shifted that they do not contaminate the 
local Ovil and Ovill absorption lines. For this reason we do not include NGC4051 in our 
final sample. We then modeled all the statistically significant AGN intrinsic absorption and 
emission features with Gaussian components. 

3.2. Local (z ~ 0) Absorption 

After fitting the continuum and intrinsic features as described above, the local O VII and 

vill Ka absorption lines are detected with > 3a and > 2a significance levels respectively. 
In 6 sources, we also detected the O VII K/3 absorption line near 18.62A. We fit these lines 
in Sherpa with narrow Gaussian features. Since with Chandra gratings (FWHM = 0.05 A 
for LEG and FWHM = 0.023 A for MEG) the lines are unresolved, we fix the line width to 

1 mA. Errors are calculated using the projection command in Sherpa, allowing the overall 
continuum normalization to vary along with all parameters for each line. For the observations 
with no detection of O VII K/3, we fixed the line centroid at 18.629 A and obtained the upper 
limits on equivalent widths (EW). The best-fit line equivalent widths (EWs) and statistical 
uncertainties are given in Table 2 and the spectra are shown in figure 1. 
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3.3. Column Density Measurement 

For optically thin gas the ionic column density depends simply on the observed equiv- 
alent width: N{ion) = 1.3 x 10 20 (-y^), where N{ion) is the ionic column density (cm -2 ), 
EW is the equivalent width (A), / is the oscillator strength of the transition, and A is in A. 
However, at the measured column densities of N(Ovil), saturation could be an important 
issue as suggested by simulations (Chen et al. 2003) and observational studies of Mrk421 
(Williams et al. 2005). Therefore to correctly convert the measured equivalent widths to 
ionic column densities, we require knowledge of the Doppler parameter b; at a fixed EW, 
column density decreases with increasing b. The low velocity resolution of Chandra gratings 
makes it unfeasible to directly measure the O VII line width. If multiple absorption lines 
from the same ion are detected, the relative equivalent widths of these lines can instead be 
used to place limits on the column density N(0 vn) and the Doppler parameter b of the 
medium. 

To use this technique, we searched for and detected Ovil K(3 line at 18.629 A in 6 
out of 8 targets and measured the upper limit for the remaining two. For O vn, the ex- 
pected ^[^-f] ratio is f^x^xx^Ka) = 0-156. Our observations indicate that most Ovil Ka 
lines are saturated (Table 2). To place quantitative constraints on N(Ovn ) and the b 
parameter we employ the technique described in details in (Williams et al. 2005). For a 
given absorption line with a measured equivalent width and known oscillator strength / 
value, the inferred column density as a function of the b parameter can be calculated 
using the relations from (Spitzer 1978). The b and N(Ovn ) can be determined for a 
range of Doppler parameters for which both Ka and K/3 transitions provide consistent 
N(Ovn ) measurements. Figure 2 shows such la contours for the measured Ovil Ka 
and K/3 transitions for NGC3783. As the figure show, the la constraints on Doppler pa- 
rameter b and Ovil column densities are 45 < b < 128 km s _1 and 16.03 <logN(Ovil )< 
16.68 cm -2 , respectively. For Mrk421 and PKS2155-304 we used the b parameter values from 
(Williams et al. 2005; Williams et al. 2007). Following the same method, we constrained the 
b parameter and Ovil column densities for all other observations. Our measured column 
density toward Mrk 421 is consistent with that in (Yao et al. 2008). The upper limit on 
the column density through the Galactic halo by Yao et al. is based on several assumption 
which, in our opinion, are faulty, but a detailed discussion is beyond the scope of this Letter. 
The Ovil column densities for our sample range from logiV(Ovil )= 15.82 to 16.50 cm~ 2 , 
with a weighted mean value of logiV(0 Vll)= 16.19 ± 0.08 cm~ 2 . These values of N(0 VII ) 
are higher than measured by other authors (Fang et al. 2006; Bregman et al. 2007) who 
assumed the lines to be unsaturated. 
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4. Results 

Column density ratios of different ions of the same element depend on the physical state 
of the medium and provide the rigorous constraints on the temperature of the absorbing 
medium. Given that we detect both O VII and O Vlll , we place a tight constraint on the 
temperature logT = 6.1 — 6.4 K, assuming the gas to be in collisional ionization equilibrium. 

With the spectral resolution of current X-ray gratings it is difficult to resolve the ab- 
sorption lines into Galactic and Local Group components. However a comparison between 
emission and absorption measurements, provides us with a great tool to analyze the prop- 
erties of the X-ray absorbers. The absorption lines measure the column density of gas 
Nh = fJ>n e R, where \i is the mean molecular weight ~ 0.8, n e is the electron density and R 
is the path-length. The emission measure, on the other hand, is sensitive to the square of 
the number density of the gas (EM = n^R, assuming a constant density plasma). Therefore 
a combination of absorption and emission measurements naturally provides constraints on 
the density and the path- length of the absorbing/emitting plasma. 

(Henley et al. 2010) and (Yoshino et al. 2009) using XMM-Newton and Suzaku data 
respectively, found that the Galactic halo temperature is fairly constant across the sky 
(1.8 — 2.4) x 10 6 K, but the halo emission measure varies by an order of magnitude (0.0005 — 
0.005 cm" 6 pc) with an average of EM = 0.0030 ± 0.0006 cm~ 6 pc, assuming solar metallic- 
ity 1 . Other measurements of Galactic Halo emission (McCammon et al. 2002; Galeazzi et al. 2007; 
Gupta et al. 2009; Hagihara et al. 2010) also reported EM ~ 0.003 cm -6 pc for solar metal- 
licity, close to the average. Thus using EM = 0.003 {^){ 8 (Ao/a h ) ) cm ~ 6 P c > we solve for the 
path length and electron density of the absorbing gas. Combining absorption and emission 
the density is: 

n e = (2.0 ±0.6 x 10- 4 )(-^)- 1 cm" 3 (1) 

JOVII 



and the path length: 



151 x 10~ 4 W 0.5 2 Z Q . 



R = (71.8 ± 30.2)( - )(-— Y(^) kpc (2) 

[Ao/Ah) JOVII A 



1 Thc X-ray spectrum of the diffuse Galactic halo emission is normally fitted with collisional ionization 
equilibrium (CIE) plasma models (e.g. APEC in XSPEC). The free parameters of this model are temperature, 
metallicity and the emission measure. In the low-resolution spectrum the emission measure and metallicity 
cannot be fit independently, so the EM is obtained for a given metallicity; smaller the metallicity, larger the 
EM. 



-7- 



where the Solar Oxygen abundance of A /A H = 8.51 x 10" 4 is from (Anders & Grevesse 1989), 
fovn is the ionization fraction of O VII and Z is the metallicity. Newer values of oxygen 
abundance are even lower (Holweger 2001; Asplund et al. 2009); making L larger. For the 
observed temperature of about > 10 6 K, it is reasonable to expect / = 0.5 (see, e.g., figure 4 
in (Mathur et al. 2003)). Cosmological simulations (Toft et al. 2002; Sommer-Larsen 2006) 
of formation and evolution of disk galaxies show that outside the galactic disk the mean 
metallicity of gas is Z = 0.2 ± 0.1 Z Q . These values of metallicities are also consistent 
with observational results for the outskirts of groups (Rasmussen et al. 2009) and clusters of 
galaxies (e.g., (Tamura et al. 2004)). Thus it is highly unlikely that in the CGM metallicity 
is as high as Z Q . As such, Z = Z Q sets a lower limit on the pathlength and an upper limit 
on the electron density. The la limits on the pathlength and density are L > 41.6 kpc and 
n e < 2.6 x 10~ 4 cm -3 . For Z = 0.3Z & , which is far more likely, the path-length becomes as 
large as L = (239 ± 100) kpc. 

We can also estimate the total baryonic mass traced by the O VII absorbers, assuming 
a homogeneous spherically symmetric system: 



M total = (2.3 ± 2.1) x lO^X 8 ; 5 ! X 7/ ! ; 4 ) 3 (^) 3 (^) 5 M ° ( 3 ) 

{A /A H ) Z JOVII 

for L = 239 kpc and f c is the covering factor, ~ 72% for our entire sample. 

For the la lower limit on the path-length (L > 139 kpc), the mass is M tota i > 6.1 x 
1O 1O M . For the Solar abundance used by (Tumlinson et al. 2011) M totai = 1.2 x 1O 1O M , 
compared to 2 x 1O 9 M found by (Tumlinson et al. 2011) in Ovi systems. Thus we find 
that the Ovil /Ovill systems probe the reservoir of gas in the CGM extending to over 
100 kpc. The mass probed by this warm-hot gas is larger that that in any other phase of 
the CGM and is comparable to the entire baryonic mass of the Galactic disk ~ 6 x 1O 1O M 
(Sommer-Larsen 2006). The baryonic fraction fl of this warm-hot gas varies from 0.09 — 0.23 
depending on the estimates of the virial mass of the Milky Way, from 1O 12 M to 2.5 x 1O 12 M 
((Anderson & Bregman 2010) and references therein), bracketing the Universal value of ft, = 
0.17. 



The oxygen mass in the CGM as traced by O VII /O VIII is: 

M oxygen = (6.8 ± 5.9) x 1O 8 (-^)M (4) 

JOVII 

and with the effective oxygen yi eld ~ 0.01 for a Salpeter IMF, the stellar mass needed 
to produce this amount of oxygen is M+ = (6.8 ± 5.9) x 1O 1O M , which is of the order of the 
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disk+bulge stars. 

4.1. Comparison with models 

(Ntormousi & Sommer-Larsen 2010) calculate all sky O VII column density distributions 
for halos of three Milky Way-like disk galaxies, resulting from cosmological high-resolution N- 
body/gasdynamical simulations. These simulations predict the mean logN(Ovil ) ranging 
from 13.46-14.55, which is significantly lower than the average column density we find. 
Our observations, however, are consistent with those of other such studies (Fang et al. 2006; 
Williams et al. 2005; Williams et al. 2007; Bregman et al. 2007). We note that in the "pure 
hot halo" models of (Ntormousi & Sommer-Larsen 2010), the Local Group contribution is 
not included; this might be a part of the reason for the discrepancy (see section 4.2). 

Our path-length and density measurements are in excellent agreement with the recent 
estimates of the distribution of hot (T ~ 10 6 K) circum-galactic gas, based on high resolu- 
tion cosmological hydrodynamical simulations of (Feldmann et al. 2012). These simulations 
predict the presence of ionized circum-galactic gas with a density of ~ 10~ 4 cm -4 out to 
100 kpc and beyond. Hydrodynamic simulations of CGM at z = by (Stinson et al. 2011) 
found that their simulated galaxies contain metal enriched warm-hot gas extending to ap- 
proximately the virial radius ( 250 kpc). While their simulations do not track the Ovil or 
Ovill producing gas, it is of interest to note that the mass in the CGM is several times 
higher than the stellar mass of their simulated galaxies, similar to what we find. They argue 
that the missing baryons exist in the CGM. 

4.2. Assumptions and Biases 

We have combined our measurements of the Ovil and Ovm column densities with 
the published values of emission measure (EM) to derive physical conditions in the absorb- 
ing/emitting plasma, viz. temperature, density, path-length and mass. The constraint on the 
temperature is robust; it depends simply on the O VII to Ovill column density ratio, under 
the assumption of collisional ionization equilibrium. For our measured density of a few times 
10~ 4 cm -3 , which is several thousand times larger than the mean density of the Universe 
(n = 1.88 x 10~ 7 cm -3 ), photoionization becomes unimportant and collisional ionization is a 
reasonable assumption. 

We have assumed that the absorbing/emitting plasma is of uniform density; this need 
not be the case. The distribution of warm-hot gas around several galaxies follows a beta- 
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model profile in which the density is high in the center and falls off with radius (see 
(Mathur et al. 2008) and references therein); gas in clusters and groups of galaxies also 
follows similar profile. In any case, if the gas is not of uniform density, the EW would 
be weighted by the denser gas and for the same observed column density the average den- 
sity would be smaller and the inferred path-length would be larger extending into the LG 
medium. We cannot rule out the possibility that the warm-hot gas we trace is from the LG 
medium instead. We note, however, that in the simulations of (Feldmann et al. 2012), the 
density is roughly constant above the Galactic disk out to about 100 kpc. 

As noted above, our sample is made of sight-lines in which both O VII and O VIII Ka 
absorption lines at z = are securely detected. Does this mean that our sample is biased 
toward high column density systems? If so, the path-length and the mass of the warm-hot 
plasma could be severely over-estimated. In figure 3 we have plotted the EW distribution 
of O VII line in all 29 sources in the parent sample. The sub-sample of 8 sources is marked 
by the red shaded region; it is clear that this sub-sample is not biased toward high EW 
systems. The high EW systems are not in the sub-sample of this paper because they do 
not have clean detections of O Vlll lines, in large part because of the contamination from 
intrinsic absorption. It is also worth noting that the upper limits on the EWs in targets with 
non-detections of O VII Ka line is higher than the average of the parent sample in all but one 
cases. All together, we argue that the average EW of our sub-sample is a fair representation 
of the large sample and that our sample is not biased toward high column density systems. 

5. Conclusions 

For reasonable values of parameters and with reasonable assumptions, the Chandra ob- 
servations of O VII and O VIII absorption lines at z = imply that there is a huge reservoir 
of ionized gas around the Milky Way. It may be in the halo of of the Milky Way or in the 
surrounding Local Group. Either way, its mass appears to be very large. 

Support for this work was provided by the National Aeronautics and Space Administra- 
tion through Chandra Award Number TM9-0010X issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of the National Aeronautics Space Administration under contract NAS8-03060. 
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Fig. 1. — The normalized flux at the location of the O VII line at 21.602 A. 
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Fig. 2. — Contours of allowed column densities N(0 VII ) and Doppler parameters b for the 
Ovn K q (black) and Ovil Kj3 (red). 
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EW 0VII (mA) 



Fig. 3. — Distribution of Ovil k« line EW for the parent sample. The solid black line 
corresponds to distribution for the 21 observations in which VII k q lines are clearly visible. 
The black shaded region marked the la upper limits on the VII k« EW for 8 observations 
with no (or less than la significance) detection. The red shaded region corresponds to 
sub-sample selected for this study (see text). 
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Table 1: Summary of the targets used in this investigation. 



Target 


1 


b 


Rcdshift 


Exposure 




deg 


(deg) 


z 


(ks) 


Mrk290 


91.48 


47.95 


0.0304 


250 


PKS2155-304 


17.73 


-52.24 


0.1160 


530 


Mrk421 


179.83 


65.03 


0.0300 


720 


Mrk509 


35.97 


-29.86 


0.0344 


460 


3C382 


61.30 


17.44 


0.0579 


120 


Ark564 


92.13 


-25.33 


0.0247 


250 


NGC3783 


287.45 


22.94 


0.0097 


905 


H2106-099 


40.26 


-34.93 


0.0265 


100 
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Table 2: The Ovil and OviII absorption line measurement. 



Target 


EW (OVU Ka) 

(ml) 


EW (Ovu k/?) 
(ml) 


EW (OVIII Ka) 

(ml) 


UVH \EW(K a )! 


b 

km/s 


log(NO vii ) 

cm~ 2 


Mrk290 


18.9 ±4.5 


5.1 ±3.7 


8.4 ±2.9 


0.27 ± 


0.21 


> 


55 


16.14 ±0.32* 


PKS2155-304 


11.6 ± 1.6 


4.2 ±1.3 


6.7 ± 1.4 


0.36 ± 


0.12 


35 


-94 


16.09 ±0.19 


Mrk421 


9.4 ±1.1 


4.6 ±0.7 


1.8 ±0.9 


0.49 ± 


0.09 


24- 


- 55 


16.22 ±0.23 


Mrk509 


23.9 ±5.0 


11.7 ± 4.1 


10.3 ±4.3 


0.49 ± 


0.20 


70 - 


- 200 


16.7 ±0.27 


3C382 


17.3 ±5.0 


7.8 ±3.0 


6.8 ±3.8 


0.45 ± 


0.22 


> 


40 


16.50 ±0.49* 


Ark564 


12.0 ± 1.9 


< 3.8 


9.5 ±4.1 






> 


20 


15.82 ±0.20* 


NGC3783 


14.4 ±2.5 


5.6 ±1.6 


4.5 ±2.9 


0.39 ± 


0.13 


50- 


- 130 


16.30 ±0.25 


H2106-099 


48.3 ± 18.0 


< 34.2 


28.8 ± 13.8 






> 


70 


16.23 ±0.16* 



*The lower limit on Ovil column densities are calculated using the curve-of-growth 
analysis. 



